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A Lifshitz scalar with the dynamical critical exponent z = 3 obtains scale-invariant, super-horizon 
field fluctuations without the need of an inflationary era. Since this mechanism is due to the special 
scaling of the Lifshitz scalar and persists in the presence of unsuppressed self-couplings, the resulting 
fluctuation spectrum can deviate from a Gaussian distribution. We study the non-Gaussian nature 
of the Lifshitz scalar's intrinsic field fluctuations, and show that primordial curvature perturbations 
sourced from such field fluctuations can have large non-Gaussianity of order /nl = O(IOO), which 
will be detected by upcoming CMB observations. We compute the bispectrum and trispectrum of 
the fluctuations, and discuss their configurations in momentum space. In particular, the bispectrum 
is found to take various shapes, including the local, equilateral, and orthogonal shapes. Intriguingly, 
all integrals in the in-in formalism can be performed analytically. 



I. INTRODUCTION 



Hofava-Lifshitz gravity ^ is attracting much attention as one of candidates for the theory of quantum gravity 
because of its power-counting renormalizability, which is realized by the Lifshitz scaling 

X bx, t — > ¥t, (1) 

with the dynamical critical exponent z > 3 in the ultraviolet (UV). There are many attempts to investigate properties 
and implications of this theory [2, 3 . 

It is natural to suppose that not only gravitational fields but also other fields exhibit the same Lifshitz scaling in the 
UV. Even if they classically have different scalings, quantum corrections should render them to have the same scaling. 
A Lifshitz scalar with z — 3 can obtain scale- invariant, super-horizon field fluctuations even without inflation [2], 
thus can source the primordial curvature perturbations through mechanisms such as the curvaton scenario [4 or the 
modulated decay [5]. It is noteworthy that this value of z is the minimal value for which gravity is power-counting 
renormalizable . 

In order to discern this production mechanism of the primordial perturbation from others, we need to investigate 
distinct features in observables such as the cosmic microwave background. In this respect, non-Gaussianity has been 
considered as one of the promising approaches to distinguish production mechanisms. For this reason, there are 
on-going efforts to detect or constrain non-Gaussian nature of the primordial perturbation [5]. Towards identification 
of the production mechanism by future observations, theoretical analyses of non-Gaussianity in various cosmological 
scenarios have been performed [TMlOj. 

In this paper, we focus on primordial non-Gaussianity from a Lifshitz scalar and calculate its bispectrum and 
trispectrum. With the dynamical critical exponent z — 3, the scaling dimension of the Lifshitz scalar is zero and, 
thus, nonlinear terms in the action are unsuppressed unless forbidden by symmetry or driven to small values by 
renormalization. It is those nonlinear terms that we expect to produce non-Gaussianity. Even when the Lifshitz 
scalar's field fluctuations are linearly transformed to the curvature perturbations (which can be realized by the 
curvaton mechanism or/and modulated decay), it turns out that the produced bispectrum can be large enough to 
be observed in future observations. We find three independent cubic terms dominant in the UV, each of which gives 
different shape dependence of the bispectrum. Roughly speaking, they correspond to local, equilateral and orthogonal 
shapes, respectively. 

The rest of this paper is organized as follows. In section ITT] we review generation of scale- invariant cosmological 



perturbations from a Lifshitz scalar. In section III we estimate the size of non-Gaussianity and see that the nonlinear 
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parameter /nl can be as large as O(IOO). In section IV we concretely show the momentum dependence of the 
bispectrum and trispectrum. Section [V] is devoted to a summary of this paper and discussions. In appendix |A] we 
derive the set of independent cubic and quartic terms dominant in the UV. 

II. SCALE-INVARIANT POWER SPECTRUM FROM LIFSHITZ SCALAR 

In this section, we review the mechanism for generation of scale-invariant cosmological perturbations from a Lifshitz 
scalar [2]. The action for a Lifshitz scalar $ in Minkowski background is 



{dt<S>f + <^d<i> + 0{<i>^) , (2) 



where 



O = (-1)"+^ . ^ A" + i^iy J'-' A"-i + • • • + siA - m\ (3) 

A = d^^didj, M and to are mass scales and s„ are dimensionless constants. Here, it is supposed that the time kinetic 
term is already canonically normalized, and thus nonlinear terms in the action indicated by 0{^^) do not include 
time derivatives. On the other hand, those nonlinear terms can include spatial derivatives. Also, the sign of the first 
term in the right hand side of (|3| is set by requiring stability in the UV. 

In the UV, the first term in O is dominant and the field $ described by the action (jsj) exhibits the Lifshitz scaling 
with 

$ ^ (4) 

We find that for z — 3, the scaling dimension of $ is zero and thus the amplitude of quantum fluctuations of $ is 
expected to be independent of the energy scale of the system of interest. This indicates that the power spectrum of 
quantum fluctuations of $ in an expanding universe should be scale- invariant. Intriguingly, the minimal value of z 
for which Hofava-Lifshitz gravity is power-counting renormalizable is also 3. Hereafter, we consider the z = 3 case. 
Now let us consider the Lifshitz scalar <I>, specialized to the case with z = 3, in a flat FRW background 

ds'^ = -dt + a{tfSijdx^dx^, (5) 

to investigate generation of cosmological perturbations. We just need to replace the volume element d'^x by d'^xa{t)^ 
and the spatial Laplacian A by A/a{t)'^ in the action (|2| with z = 3. We expand the scalar field $ around a 
homogeneous v.e.v. $o as $ = $o + Throughout this paper we consider the UV regime in which the Hubble 
expansion rate H is much higher than mass scales in the scalar field action. In this regime, the Hubble friction is so 
strong that the time dependence of the background $o is unimportant. For this reason, hereafter, we treat $o as a 
constant. The action for the perturbation (f> is then written as 



dtd^x a{tf [(9t0)2 -I- (j)0<j) + 0(0^)] , (6) 



where 



= ^ A^ - A^-h-^A-TO^ (7) 

M and m are mass scales and s and are dimensionless constants. In the UV, the quadratic action for is simply 

S2^\j dtd^x a{tr I + -^^A^^l . (8) 

As discussed after Q, the scaling dimension of $ and thus is zero, 

cj, -> &>, (9) 

and its power-spectrum should be scale- invariant. Since (j) is scale-invariant and there is only one scale M in the UV 
quadratic action ([s]), we expect that the power-spectrum should be roughly 

- M^. (10) 
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FIG. 1: A schematic picture of the evolution of cosmological perturbations for a power-law expansion a oct^ with 1/3 < p < 1. 
The physical wavelength exits the sound horizon ~ {M^ H)^^^^ in the UV and re-enters the sound horizon ~ in the IR. We 
suppose that the scale-invariant perturbations of the Lifshitz scalar are converted into those of radiation by curvaton mechanism 
or/and modulated decay when physical wavelengths of interest are outside the sound horizon. Thus, strictly speaking, the sound 
horizon in the UV is for the Lifshitz scalar and that in the IR is for radiation. 



Now let us calculate the power spectrum concretely. By solving the Heisenberg equation obtained from the quadratic 
action (IS]), operator can be expanded as 



and 



where 



M f dt' \ 



(11) 
(12) 
(13) 



The mode function Uf,{t) is chosen so that its asymptotic behavior in the Minkowski limit (a(t) — > const.) is the same 
as the positive- frequency mode function in Minkowski background. The vacuum state |0) is defined as usual by 

flfc |0) = 0. (14) 

The mode function Uk{t) approaches a constant value in the a{t) — > oo limit if the integral J*°° dt/a{t)^ converges, 
where ioo is the time corresponding to a{t) — )> cxi. The power-law expansion a{t) cx with p > 1/3 satisfies this 
condition. Under this condition, when the physical wavelength a{t)/k becomes as long as the size of the sound horizon 
~ {M'^H)^^/^ , the mode function Uk{t) stops oscillating and freezes out. (See in Fig. ([I]).) Note that the physical 
wavelength at sound horizon exit is super-horizon size in the UV, i.e. when H ^ M . 
The commutation relations of operators are defined in the usual manner as 



[0(a;), 7r(y)] = i5'^{x - y), [0(a;), 0(y)] = [7r(a;), 7r(y)] = 0, 



(15) 
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where 



7t{x) = 



SS2 



^a{tfdt<l>{x), 



5dt(t){x) 

and hence the operators and a\, satisfy 

K, al,] = {2TTf5^{k - k'), [ofc, ay] = [a^.a^. 
Power spectrum is defined as 

(O|0fc0fe'|O> = {2nf5\k + k')P^, 

so that 



2^ 

1^ 



(010^10) 



dk. 



Substituting eqs.(f2) and (f3) into cq.(f8), we obtain 



(16) 



(17) 



(18) 



(19) 



(20) 



In this paper, we suppose that the primordial fluctuations of the Lifshitz scalar field (j) are almost linearly transformed 
to the curvature perturbations and that large non-Gaussianity is not created in the transformation itself. In scenarios 
like the curvaton mechanism and the modulated decay fS] , this assumption can be realized naturally, while violation 
of the assumption is also possible in these scenarios. With the assumption, the curvature perturbation ^ is related to 
the field fluctuation 6 as 



C = -+0(1) X 



0(1) X 



(21) 



where /i is an energy scale. Therefore, the power spectrum of the curvature perturbation takes over the scale-invariance 
as 



M2 



(27r)2/i2- 



The COBE normahzation |6] sets T'l''^ ~ 4.9 x 10"^ and thus 



M 



3.0 X IQ-^ 



(22) 



(23) 



III. ORDER ESTIMATE FOR /nl 



Let us move on to discussions about non-Gaussianity. In the previous section we have seen that a free Lifshitz 
scalar with the dynamical critical exponent z = 3 can produce scale-invariant perturbations. The essential reason for 
this is that the scaling dimension of the scalar is zero at high energy. This implies that at high energy, some nonlinear 
operators could be as important as leading operators in the quadratic action and that large non-Gaussianity could 
be generated by those nonlinear operators. On the other hand, because of power-counting renormalizability, those 
nonlinear operators do not completely dominate over the quadratic terms and thus do not spoil the analysis of 
scaling dimensions and the scale-invariance of the power-spectrum unless the theory gets really strongly coupled. In 
this section, we shall present order estimates for the bispectrum of curvature perturbations and the corresponding 
nonlinear parameter /nl, deferring detailed calculations until the next section. 



As presented in (21 1, throughout this paper we assume that perturbations of the Lifshitz scalar are almost linearly 



transformed to curvature perturbations. In our calculation of the bispectrum and trispectrum of curvature perturba- 



tions, we thus ignore nonlinear terms in the right hand side of (21) and take into account the linear term only. In 
particular, 

(CO ^ M"' (00) , (CCC) ^ M"' (000) ■ (24) 
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This treatment is justified by the fact that perturbations of the Lifshitz scalar have large non-Gaussianity. Actually, 
the bispectrum and trispectrum would obtain corrections from nonlinear terms in (21 1, but those corrections would 



be smaller than large contributions from the Lifshitz scalar's intrinsic non-Gaussianity. In terms of the nonlinear 
parameter /nl, corrections due to nonlinear terms in ( plj ) are typically 0(1) but we shall see below that the scalar 
field's intrinsic non-Gaussianity can lead to /nl = O(IOO). 

We shall adopt the so called in-in formalism [7l [11] to calculate the bispectrum and trispectrum of the Lifshitz 
scalar. The leading contribution to the bispectrum is given by the following formula (see the next section for details) 



(25) 



where represents cubic terms in the interaction Hamiltonian. Dominant terms in J dtH^ are marginal ones, i.e. 
those terms whose scaling dimensions are zero. Actually, as shown in Appendix |A] there are three (and only three) 
independent marginal cubic operators in the action in the UV: 



dtd^a 



M^a{t)' 



(26) 



where ai are dimensionless parameters. (The fist term can be forbidden by the shift symmetry if one likes.) Evidently, 
validity of perturbative expansion (in the in-in formalism) requires be smaller than unity. The corresponding cubic 
operators in the interaction Hamiltonian are 



1 



d^'x {ai02A30 + a2(A20)(a,0)2 + a^iAcj,)'} 



(27) 



Each of these dominant cubic terms includes six spatial derivatives and gives zero scaling dimension to J dtH^. 
Combining this with the fact that the scaling dimension of (j) is zero, we conclude that the bispectrum of (j) given by 
(25) should be scale-independent, and thus 



where a stands for the most dominant one among ai {i = 1,2,3). 
Roughly speaking, the non-linear parameter /nl is defined so that 



/i 



NL 



(CCC) 
(CO'' 



Thus, combining this with (24 1, (10 1 and (28), we obtain 

'M' 



/nl ^ a 



3 X 10-^a. 



(28) 



(29) 



(30) 



Here, we have used the COBE normalization (23). In the next section, we obtain a more precise expression for / 



NL 



in terms of ai. 

As already stated, validity of perturbative expansion requires that the dimensionless parameters be smaller than 



unity. We find from the order estimate (30) that /nl can be large, e.g. as large as O(IOO), even if ai are reasonably 
small. 



IV. SHAPES OF NON-GAUSSIANITIES 



We have seen that the Lifshitz scalar's non-Gaussian intrinsic fluctuations can leave observable non-Gaussianities in 
the sky. In this section, we compute the bispectrum and trispectrum of Lifshitz scalar fluctuations and discuss their 
shapes. Since a Lifshitz scalar with z = 3 can have all possible self-coupling terms containing six spatial derivatives 
(cf. (26), (|53[)), the resulting correlation functions can take various conflgurations in momentum space. We will 



see especially that the generated bispectrum includes local, equilateral, and orthogonal shapes. We also discuss non- 
Gaussianity of the primordial curvature perturbations sourced by the Lifshitz scalar's field fluctuations. By comparing 
with the latest observational constraints, we obtain bounds on the Lifshitz scalar's self-coupling strengths. 
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A. Bispectrum 



We make use of the prescription of the in-in formahsm (see e.g. [71 [TT]) for calculating expectation values of a 
product Q{t) of field operators at time t, 



{Q{t)) = 



Texpii / Hi{t')dt' 



Texp 



'I I Hi{t')dt' 

'in 



(31) 



where T denotes time-ordering, T is anti-time-ordering, [t) is the product Q{t) in the interaction picture, and Hj{t) 
is the interaction part of the Hamiltonian in the interaction picture. We take the time to to be at very early times when 
the fluctuation wavelengths are well inside the sound horizon. In the present case, introducing the time dr = dt/a^, 
then the integral can be taken in terms of t from — oo to some time after the sound horizon exit.^ 

In order to obtain the bispectrum of the Lifshitz scalar fluctuations, we expand (31 ) in terms of Hj and compute 



terms which are first order in the 3-point interaction Hamiltonian (27). We then find 



(0fc,(<)0fe,(<)0fe3(t)) =i / dt'{[H3{t'),(l)kAt)<l>kAt)^ks{t)] 



dt' 



(27r)3M3 (5(3)(fci + fca + fc3)/(fci, fc2) 
22 {kik2k3)^{kl + k^ + kl) ' 



+ c.c. 



(32) 
(33) 
(34) 

(35) 



where we have defined 



fiki, kj) ^ 2ai(fcf + k] + fcf_,.) + a2(fcf -f- k) + k% - km - kfk^ 
with kij = \ki + kj\. 



kjkt - kfkl - k%t - ^kl) + 6a3kfk]k^j, (36) 



In order to study its configuration in momentum space, let us express the bispectrum ( 35 1 as 
(0fejt)0fc.(i)0fe3(^)) - (27r)3Af3^(3)(fci-f fe2 + fc3)m,fc2,fc3), 



(37) 

and plot (fci fc2fc3)^-F(fci, fc2, fcs) as a function of X2 = ^2/^1 and 3:3 = k^/ki. Contributions from the ai, a2, and 
terms in ( 36 ) are plotted respectively in Figures [2] [3j and |4] We have assumed < 0:3 < 0:2 < 1 to avoid showing the 
same configuration twice, and X2 > 1 — further follows from the triangular inequality. 




Note that the lower Umit of integration in ^3l\ is shifted to — oo(l — ie) on the right side of , and to — oo(l + it) on the left side, so 
that the oscillating exponents become exponentially decreasing. This prescription corresponds to picking up the vacuum state. 



1.0 



FIG. 4: The shape of the 0:3 contribution to 
{kik2k'if F , where Q3 — —1. 

Here, one should note that only the contribution from the ai term blows up in the "squeezed" triangle limit, i.e. 
A;3 <C ^2 ~ fci, giving a shape similar to that of the local form [12 

i^local(fcl,fc2,fc3) = 2(p^ + P^ + p'p). (38) 



This can be understood as follows: Among the three interaction terms in (26), only the ai term breaks the shift 
symmetry <j) ^ <j) -\- const. Such a shift symmetry makes the theory indifferent to infinitely large scale fluctuations, 
and suppresses interactions with long wavelength modes. This is why {kik2k^y'F from the a2 and contributions 
disappear in the squeezed limit. 

Moreover, the shapes from the a2 and terms respectively look similar to that of the orthogonal [13 and equi- 
lateral |14j forms which are given by 

/ 3 3 3 8 3 \ 

i^orthog.(fci, fc2, fcs) = 6 -pp - pp - pp - ppp + ppp + (Sperm.) , (39) 

\^ /\'-|^n,2 ni-ynj'^ i\i'2'^'^ 3 l 2 3 / 

J^cquii.(fci,A:2,fc3) = 6(~Ti:3-p^-p^ - p^ + + (Sperm.)) . (40) 

The permutations act only on the last terms in the parentheses. 



To give a bit more quantitative discussion on shapes, let us express (3S) in terms of the above forms. Following [131 



I1S| . we introduce a scalar product between two distributions Fi and F2, 

»1 rl 



Fi-i^2=> OC / dx2 dx^X2X^Fi{l,X2,X'i)F2{l,X2,X'i), (41) 

-ffci ^fc2 -ffea "'0 JI-X2 



where summation is taken over all fcj's which form a triangle, and o< k denotes the power spectrum. Then one 



can "expand" F(fci,/c2,fc3) of (37) in terms of (38), (39), (40) and obtain a template function 



^tcmplate(fcl, ^2, ^3) = (^^^^ Fiocsxiki , fc2, ^3) + C^l*""*^' -^orthog. (^1 , ^2 , ^3 ) + C^T'' ■f'oquil. (^1 , ^2 , ^s), (42) 

where 

= J'^V = -0-125ai, (43) 

-''^local ' -^local 

c°^'^°^- = ^ ' = 0.226ai + 0.0186^2 - 0.00334a3, (44) 

-^orthog. ' -^orthog. 

^equii. ^ J^'-^'^q^"- = -0.223ai + 0.0280a2 - 0.0876a3. (4S) 

-^equil. ' -^equil. 
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a2 and as are absent in (43 1 since the denominator Fiocai • -Fiocai blows up.^ Here, since (38), (39 1, and (40) do not 



form a complete basis set, one should consider (|42| as an indicator roughly telling which combination of the 3-point 
interaction terms in ( 26 ) gives which bispectrum shape. 

Using the above template function, one can estimate the non-Gaussianity parameter /nl of the primordial curvature 
perturbations. As discussed around (24), we assume that the Lifshitz scalar fluctuations are linearly converted to the 
curvature perturbations, i.e. 



c 



(46) 



with some mass scale fi. For the bispectrum of Bardeen's curvature perturbations ^I^ (which is related to the primordial 

(47) 



curvature perturbations by C = f'J' in the matter-dominated era 



(*fc,*fc.*fc3) = (2^)3<5(3'(/Ci+fc2+fe3)i^*(fcl,fc2,fc3), 

we define its non-Gaussianity parameter /nl as 



F^{k,k,k) ^ hi^Al 
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Here, A.p is the amplitude of the power spectrum 



(*fci*fc.> = (2vr)^'5(^^(fei+fc2)^: 



which in the present case is 



A* = — 



9 fM 



50 V M 



1.7 X 10" 



(48) 



(49) 



(50) 



The value on the far right hand side is set by the COBE normalization (23). Then, by substituting the template 
function (42) into the Lifshitz scalar fluctuations (37 1, we obtain 



/nl — ~3~^^NL — 2.2 X lO^c^L, 



(51) 



where i = local, orthog., equil. We can now use this relation to set constraints on the self-coupling strengths, by 
computing the necessary conditions for ( 43 1 , ( 44 ) , and ( 45 ) to satisfy the WMAP 7yr 6J constraints on the non- 



Gaussianity parameters: -10 < f^l""^ < 74, -410 < /^°^- < 6, and -214 < /^f < 266 (each at 95% CL), 
respectively. Thus we arrive at 



0.03 < ai < 0.004, 



-1 < a2 < 0.4, 



-0.5 < ag < 0.3. 



(52) 



We expect that more rigorous treatments (such as using observational data to constrain combinations of different 
shapes, instead of individual ones) would not change the results significantly.^ 

The bounds in (52) show that observational constraints require \a2\, jaaj = 0(10^^). Note that such smallness of 
the coupling strengths is required anyway for validity of the perturbative expansions in terms of Hj, which we have 
carried out in computing the bispectrum. In other words, as a2 or become of order unity and our procedure for 
computing the bispectrum breaks down, the curvature perturbations saturate the current observational limit for the 
orthogonal/equilateral form bispectra. The most stringent bounds are obtained for ai, i.e. |ai| = O(10~^ — 10~^) 
where the required level of tuning depends on its sign. This is mainly due to the rather tight constraints on the local 
form bispectrum it produces. In order to suppress ai, some sort of symmetry may be required in the theory. For 
example, as we have stated below (38), a shift symmetry for forbids such self-coupling terms producing local- type 
bispectra. 

Here we have focused on non-Gaussianity sourced by the intrinsic fluctuations of the Lifshitz scalar, but we should 
remark that further non-Gaussianities in the resulting curvature perturbations can be generated through non-linear 
conversion processes. 



^ It is possible to avoid such divergences by introducing cutoffs in the integration, e.g., X2,X3 > 0.001 which roughly corresponds to the 

current observable limit in the CMB. 
^ We also note that ( |52[ | are conservative bounds, since they are necessary conditions for satisfying observational constraints. One way to 

obtain more strict bounds is by imposing observational constraints under an ad hoc assumption that two out of the three a^'s are zero. 

However, such procedure does not change the results significantly, except for the strong upper limit on f^^^^°^' tightening the upper 

bound on Q2 and lower bound on «3 roughly by an order of magnitude. 
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B. Trispectrum 



The trispectrum can be obtained in a similar manner to the previous subsection. We compute contributions from 
the scalar-exchange diagram (Figure [5]), and from the contact-interaction diagram (Figure [6]). The 4-point interaction 
terms of the Lifshitz scalar action are given in Appendix [X] as 



54 = / dtd-^x 



+/340'(Aa,0)2 + Pr^cj.^d^d.dkcpf -f /36(5,9,5fc0)(a,0)(9,0)(afc0)}, (53) 
where (3i are dimensionless parameters. Therefore the 4-point interaction Hamiltonian becomes 

+/3402(Aa,0)2 + k^\d,d,dk<i>f + i3^{d,d,dkc^m(^){d,(j)){du(^)]. (54) 

Even if the 3-point interactions are suppressed, i.e. \ai\ <^ 1, the 4-point interaction terms can produce a large 
trispectrum through the contact-interaction diagram. 






FIG. 5: The scalar-exchange diagram. 



FIG. 6: The contact-interaction diagram. 



The contribution from the scalar-exchange diagram are obtained as follows (Note that we only compute contributions 
from connected diagrams.), 



{4>kAt)<t>kAt)<t>k.{t)^kM)sc = \~ I dt' I dt''{[H3{n,[H3{t'),<i>kAt)<i>k,mksmkAm) 

\ Jto Jto 



connected 



(27r)3M4 (5(3) (fci + k2+k3 + ki){kl + 



kl + kl^)f{ki,k2)f{ks,k^) 



(55) 



(56) 



23 {kik2kzkiki2f{kl + kl + kl2){kl + fc| + hil^){kl + kl + kl + kl) 
+ (2 terms with (1,2)(3,4) (1,3)(2,4), (l,4)(2,3)y (57) 



where f(ki,kj) is defined in (36). One can also compute the contribution from the contact-interaction diagram, 



{^kAt)<pkAt)<i>ks{t)<i>kAt))c.i.^ r / dt'{[Hi{t'),ci)kAtHk,mksmkAm 



ta 



connected 



{2n fM^ S^^\ki + k2 + ks + ^4)^(^1, fca, fcg, k^) 



23 



{kik2k3kiY{kl + kl^kl + kf) 



(58) 
(59) 
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where we have defined 



r(fci,fc2,fc3,fc4) = 

6/3i {kf + kl + kl + 

~|~ '^ijS^ ^^1^2 ~^ ^1^2 ^1^3 ^1^3 ~^ ^1^4 k-^k^ ~j~ ^3 ~f~ ^2^3 ~^ ^2^4 ^2^4 k^^k^ ~\- ^3^4) 

~|~ 6^3 {^1 ^2 ^3 ^1 ^2 ^4 k^k^k^ -|~ k')k':ik^ 



4^4 {kikl{k^ ■ fca) + fc?fc|(fci • fcg) + fc?fc|(fcl • ^4) + fc^fc2(fc2 • fcg 



2 1,2/ 



+ 4/35 ((fci • fca)' + (fei • fca)^ + (fci • fc4)-'' + (fc2 • fes)' + (fc2 • 

+ 6/36 ((fcl • fc2)(fel • fc3)(fcl • fc4) + (fe2 • fcl)(fe2 • fc3)(fc2 ' ^4) 
+ (fca • fcl)(fc3 • fc2)(fc3 • ^4) + (^4 • fcl)(fe4 • fe2)(fe4 

Expressing the trispcctra as 



+ klkl{k2 ■ k^) 
(fe3 • fc4)') 



klkl{k^-k^)) 



fea)). 



{M)Mt)4>ks{t)4>kM) = (2^)'MV4)(fci +fc2 +fe3 +fc4)T(fei,fc2,fc3,fc4), 



(60) 



(61) 



let us plot (fcifc2fc3A:4)^/^T(fci, ^2, fc3, fe4) in the "equilateral" limit: ki = k2 = k^ = k4^. In this limit, the shape of the 
tetrahedron formed by fcj's depends on two independent variables, e.g., C2 = fci • k2/kik2 and C3 = fci • k^/kik^. 
Contributions to 7^ o. Ist] ) from the a\, 0:3, aiQ;2, aiQ!3, and a2a3 terms are plotted respectively in Figures [vj [sj 
[9j[l0j[TT] and 12 Contributions to TcA. (59) from the /3i, terms are plotted in Figures 13 Tlj and 15 We have 

omitted contributions from the (32, Ps, Pi terms since in the equilateral limit they become flat and are equivalent to 
that from /3i up to overall constant factors. The condition C2 + C3 < is required for fcj's to close, and we further 
assume C2 < C3 in order to avoid showing the same configuration twice. 




FIG. 9: The shape of the Ofl contribution to 
{kik2k3k4)^^^Ts.B., where a| = 1. 



FIG. 10: The shape of the 0:10:2 contribution to 
(^1/22^3^4)^'' ''7^. c, where 0102 = — 1- 



FIG. 11: The shape of the aias contribution to 
{k-ik2k3kif^^Ts.c., where Q1Q3 — 1. 



FIG. 12: The shape of the 020^3 contribution to 
{kik2k2,k4f^*Ts.c., where a2az = —1. 




FIG. 15: The shape of the j3& contribution to 
(^1^2^3/04)®^"*?^.!., where = —1. 
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V. CONCLUSIONS 



In this paper, we have studied non-Gaussianity in the intrinsic fluctuations of a Lifshitz scalar which follows an 
anisotropic scaling with z — S. Our work is based on [2J, which pointed out that its special dispersion relation in the 
UV can lead to generation of super-horizon field perturbations. Since the scaling dimension of a Lifshitz scalar with 
z = 3 is zero, the resulting field perturbations become scale-invariant whether or not the scalar's self-couplings are 
small. This leads to our main point that curvature perturbations generated from such field fluctuations necessarily 
leave large non-Gaussianity in the sky, unless the field's self-couplings are forbidden by some symmetry, or the field 
exhibits some sort of asymptotic freedom. This is to be contrasted with perturbations generated through cosmic 
inflation, where largely non-Gaussian intrinsic fluctuations are in most cases incompatible with scale-invariance. 

The Lifshitz scalar's self-coupling terms containing spatial derivatives produce non-Gaussianities with various con- 
figurations in momentum space. In particular, the bispectrum of the field fluctuations includes shapes which are 
similar to that of the local, equilateral, and orthogonal forms. (However, we emphasize that the local, equilateral, 
and orthogonal shapes do not form a complete basis set for the bispectrum obtained in this paper. We also note 
that the results of the effective field theory approach in [16] do not apply to our case, where Lorentz symmetry is 
explicitly broken, and non-Gaussianity is sourced by marginal terms in the action.) Upon computing the correlation 
functions, we have carried out expansions in terms of the interaction Hamiltonian. Within the domain of applicability 
of such perturbative expansion, i.e. the self-couplings less than unity, we have seen that the Lifshitz scalar's field 
fiuctuations can lead to significant non-Gaussianity in the primordial curvature perturbations. In particular, when 
curvature perturbations are sourced linearly from the field fluctuations as in (461, their bispectrum saturates the 
current observational limit for the orthogonal and equilateral forms, as the self-couplings a2 and ag, in (26) approach 
unity. Since naively there is no reason for such self-couplings to be suppressed, we can expect large non-Gaussianity to 
be produced from Lifshitz scalar fluctuations, which may be detected by upcoming CMB observations. On the other 
hand, for the local-type bispectrum, observational constraints require ai to be as small as O(10~^— 10"'^) (the level of 
tuning depends on ai's sign). However, as we have remarked, such self-couplings sourcing local-type non-Gaussianity 
can be forbidden by a shift symmetry. 

The field fluctuations generated in the mechanism of [5] obtain a scale-invariant spectrum. However, when one 
takes into account the renormalization-group fiow of the parameters of the theory (e.g. M in ([s])), the spectrum may 
become tilted. A time-dependent background value $0(0 may also give rise to similar effects. How strong the tilt 
becomes, as well as the scale-dependence of the non-Gaussianity, remains to be understood. While in this paper we 
have studied fluctuations of scalar flelds, the scalar graviton which can show up in Hofava-Lifshitz gravity may also 
obtain fluctuations in a similar manner. It would be interesting to investigate the possibility that such scalar graviton 
generates the primordial curvature perturbations. (Ref. |17j works in this direction. See jlSj for some issues related 
to the scalar graviton, including non-perturbative continuity of the limit in which general relativity is supposed to be 
recovered. See also fl9' for a recent attempt to eliminate the scalar graviton from the theory.) Furthermore, when 
considering cosmic inflation in Hofava-Lifshitz gravity, due to the field fiuctuations freezing-out at the time of sound 
horizon (A-PH)^^/^ exit, the well-known relations in slow-roll inflation between various cosmological observables and 
the slow-roll parameters are expected to be modified. Aspects of cosmic inflation in Hofava-Lifshitz gravity are also 
worthy of study in details. 
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Appendix A: Nonlinear terms in the UV action 



In this appendix we construct nonlinear, marginal terms in the action in the UV, specializing to the z = 3 case. 
As mentioned after (jsj), those nonlinear terms do not include time derivatives but can include spatial derivatives. 
We demand that each term in the action can include only up to six spatial derivatives. This treatment is justified 
since with z = 3, terms with more than six spatial derivatives would be non-renormalizable and thus would not be 
generated by quantum corrections. In this case the most important terms in the UV are marginal ones, i.e. those 
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with six spatial derivatives. In the following we write down all independent combinations of cubic and quartic terms 
which are marginal in the UV ^. 



1. Cubic terms 



In general we can write down fourteen cubic terms with six spatial derivatives as follows, 



A = (Aa,0) (Aa,(/.) 0, = (d^djdkcj)) {d,djdk4>) 0, A^ 

= (Aa,0) (A0) , Alo = (Ac',^) (9,0) 
^12 = (A0) (A0) (A0) , Ai3 = (A0) (9,9,0) (9,9,0) 



A3 = (A20) (A0) , = (A9,9,0) (9,9,0) 0, 

A20) (9,0) (9,0) , As = (A9,9,0) (9,0) (9,0) , (Al) 
All = (9i9,9fc0) (9i9,-0) (9fe0) , 
Ai4 = (9i9,0) {djdk4>) (dkdicj)) . 



Some linear combinations of these terms turn out to be total derivatives. Thus, not all of these fourteen terms lead 
to independent terms in the action, provided that surface terms are dropped. Apparently, there are thirteen total 
derivatives with six spatial derivatives as follows. 



9, [(A29,0)02] = A1 + 2A2, 
9, [(A9,9,0) (9,0) 0] = A2 + A4 + As, 
9, [(A9,0) (9,9,0) 0] - A4 + As + Aio, 
9, [(A9,0) (9,0) (9,0)] = A7 + 2Aio, 
9, [(9,9,9fe0) (9,0) (9fe0)] = As + 2Aii, 
9, [(A0) (9,9,0) (9,0)] = Ag + Aio + A13, 
9i [(9,9fe0) (9,-9fe0) (9,0)] = 2An + A13. 



9, [(A20) (9,0)0] -A2 
9, [(A9,0) (A0) 0] = A3 - 
9, [(9i9,-9fe0) (9,-9fc0) 0] = 
9, [(A9,0) (9,0) (9,0)] 



f A3 + A7, 

I- As + Ag, 

= A4 + Ag + All, 

A8 + A9 + A10, 



9,[(A0) (A0) (9,0)] = 2Ag + Ai2, 
9, [(9,9,0) (9,9fc0) (9fc0)] = Aio + An + A 



(A2) 



14, 



Actually, two of the combinations in ( A2 ) are not independent and ( A2 ) include only eleven independent combinations. 
As a result, we can represent arbitrary marginal cubic terms in the action as linear combinations of three among 
fourteen in ( Al ). 



Concretely, in this paper we shall use the following three terms. 



5*3 = / dtd 



1 

M5a(t)3 



{ai02A30 + a2(A20)(9,0)2 + a3(A0)3} 



(A3) 



An advantage of this representation is that contributions from these three terms roughly correspond to local, equilateral 
and orthogonal shape of bispectrum, respectively. 

If we impose the shift symmetry for then the number of independent interaction terms is reduced. With the shift 
symmetry, we have to consider only the last eight terms in (All and only the last seven combinations in (A2). One 
of the last seven combinations in (A2| is not independent. As a result, the marginal cubic terms in the action with 
the shift symmetry can be written as (A3) with ai = 0. 



2. Quartic terms 



Similarly, we can write down twenty quartic terms with six spatial derivatives as follows. 



Bl E 


e(A30) 03, 




B2E 


E (A29,0) (9,0)02, 


S3 EE 


(A20) (A0)02, 






E (A9,9,0) (9,9,0) 02, 




^5 = 


E (A20) (9,0) (9,0) 0, 




(A9,9,0) (9,0) (9,^ 


/')'/' 


B7 = 


e(A9,0) (A9,0)02, 




i?8 = 


E (9,9,9fc0) [didjducj)) 0^, 


Bg = 


(A9,0) (A0) (9,0) 


^, 


Bio 


= (A9,0) (9,9,0) (9,0 




Bii 


= (didjdkfp) {didjcj)) {dkcf) 0, 


B12 E 


E (A0) (A0) (A0) 




Bi3 


= (A0) (9,9,0) (9,9,0 




Bi4 


= {didjcj}) (9,9fc0) (9fe9,0) 0, 


-Bis = 


E (A9,0) (9,0) (9,0; 


\dj 


BlQ 


= (9i9,-9fc0) (9i0) (9,-0 


) (9^0) 


, Bl7 


= (A0) (A0) (9,0) (9,0) , 


Bis = 


E (A0) (9,9,0) (9,0; 


(9, 


Big 


EE (9i9,-0) (9,9,0) (9fe0 


) (9^0) 


, S20 


= (9,9,0) (9,9fe0) (9,0) (9fe0) 









In [20) . non-Gaussianity of cosmological perturbations in Hoi'ava-Lifshitz gravity is discussed. However, in that paper, non- 
renormalizable interaction terms are considered as the main source of non-Gaussianity. In the present paper, on the other hand, 
all terms in the action are power-counting renormalizable and we consider marginal ones as the main source of non-Gaussianity. 
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Apparently, there are sixteen total derivatives as follows, 



a31 _ 



= Bi+ 352, 



= B2 + B4 + 2i?e, 

= i34 + i?7 + 2i?io, 

(A9,0) (9,0) (9,0) 0] = B5 + 2Bio + Bi5, 
(9,9,9fc0) (9,0) (9fe0) 0] = -Be + 2Bn + B,e, 



(A0) (9,9,0) (9,0) 0] = + Bio + -B13 + B18, 5, 
(9,9fc0) (9,9fe0) (9,0) 0] = 2Bn + B13 + Big, 9, 
(9,9fc0) (9,0) (9,0) (9fe0)] = B16 + Bis + 2B20, 9, 



(A20) (9,0)02] =B2 + S3 + 2B5, 
(A9,0) (A0)02] =B3 + B7 + 2B9, 
(9,9,9fe0) (9,9fc0) 02] = B4 + + 2Bii, 
(A9,0) (9,0) (9,0) 0] = Be, + Bg + Bio + Si5, 
(A0) (A0) (9,0) 0] = 2B9 + B12 + Bi7, 
(9,9,0) (9,9fc0) (9fc0) 0] = Bio + -Bii + Bu + B20, 
(A0) (9,0) (9,0) (9,0)] = Bi5 + Bi7 + 2Bi8, 
(9,9,0) (9,0) (9fe0) (9fc0)] = Bi5 + Big + 2B20. 



(A5) 



Two of them are not independent. Thus, marginal quartic terms in the action can be written as a linear combination 
of six independent terms. In this paper, we use the following terms. 



*S'4 



dt(fx 



1 



M6a(i)= 



{M^A^ + /320'(A0)(A20) + /330(A0)3 

+/340'(A9,0)2 + /3502(9,9,9fe0)2 + /36(9,9,9fc0)(9,0)(9,0)(9fc0)}. (A6) 



With the shift symmetry, we have to consider only the last six terms in ( A4 ) and only the last three combinations 
in (A5|. Therefore, the number of the independent terms is three. In the notation (A6), the shift symmetry imposes 

3 



the constraints as 



/3i 



1 

12' 



/32 =/34 + ^/35, 



-/34-:^/?5, 



(A7) 



and then the action can be transformed into the following form. 



^4 



dtd-^x ■ 



1 



M6a(i)3 

This manifestly has the shift symmetry 



{(/34 + 2/35) (A0)' (9,0)' + /35 (9,9,0)' (9fe0)' + /?6(9,9,9fc0)(9,0)(9,0)(9fc0)}. (A8) 
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